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The combined effects of hole doping and magnetic dilu- 
tion on a lamellar Heisenberg antiferromagnet are studied in 
the framework of the frustration model. Magnetic vacancies 
are argued to remove some of the frustrating bonds generated 
by the holes, thus explaining the increase in the tempera- 
ture and concentration ranges exhibiting three dimensional 
long range order. The dependence of the Neel temperature 
on both hole and vacancy concentrations is derived quantita- 
tively from earlier renormalization group calculations for the 
non-dilute case, and the results reproduce experimental data 
with no new adjustable parameters. 
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Since the discovery of high-T c superconductors much 
effort was invested in the investigation of the effect of 
dopants on the magnetic properties of the parent com- 
pounds La 2 Cu0 4 (LCO) and YBa 2 Cu 3 6 (YBCO). It is 
now well established that even a very small dopant con- 
centration, which introduces a concentration x of holes 
into the CuC>2 planes, strongly reduces the Neel temper- 
ature, Tn. In LCO, doped with strontium or with excess 
oxygen, the antiferromagnetic long range order (AFLRO) 
disappears at a hole concentration x c ~ 2% [|l], while in 
YBCO x c w 3.5% In contrast, the effect of Cu 

dilution by nonmagnetic Zn is much weaker. Like in per- 
colation, the AFLRO persists at Zn concentration, z, as 
large as 25 % {§. 

Recently, Hiicker and coworkers B studied the phase 
diagram of La2_ a; Sr. z ,Cui_2Zn z 04, and found surprising 
results: It appears that the vacancies introduced by Zn 
doping weaken the destructive effect of holes (introduced 
by the Sr) on the AFLRO. E. g., in a sample with 
z = 15%, the critical concentration x c of the holes is 
approximately 3%, i.e. larger than in vacancy free LCO. 
Also, at x = 0.017 the Neel temperature has a maxi- 
mum as function of z, implying a reentrant transition! 
To explain these phenomena, Hiicker et al. measured the 
variable range hopping conductivity in their samples (at 
temperatures lower than 150 K all samples were insula- 
tors), and showed that Zn doping lowers the localization 
radius of the holes. Their qualitative conclusion was that 
as the holes become more "mobile", their influence on 
T/v increases. However, so far there has been no quanti- 
tative understanding of the combined dependence of T/v 
on both x and z. 



In this paper we present a quantitative calculation, 
which reproduces all the surprising features of the func- 
tion T/v(x, z). Our theory extends an earlier calculation 
|| , which treated the effects of quenched hole doping on 
the AFLRO in Sr doped LCO, i. e. calculated T N (x,0). 
The same parameters were then used to reproduce the ob- 
served T/v(x) for the bi-layer Ca doped YBCO M. Here 
we reproduce the full function Tn(x, z), with practically 
no additional adjustable parameters. 

Our theory is based on the frustration model [Q , which 
argues that when a hole is localized on a Cu-O-Cu bond 
H||, it effectively turns the interaction between the Cu 
spins strongly ferromagnetic, causing a canting of the 
surrounding Cu moments with an angle which decays 
with the distance r as 1/r. The frustrating bond thus 
acts like a magnetic dipole |7|,|l0|. As argued in Ref. ||, 
similar dipolar effects also arise when the hole is local- 
ized over more than one bond. The frustration model 
also predicted a magnetic spin glass phase for x > x c 
JtJ, as recently confirmed in detail in doped LCO and 
YBCO J§[y]-|l3). Furthermore, the model successfully 
reproduced the local field distributions observed in NQR 
experiments [T^ ] . In earlier work, Glazman and Ioselevich 
]l5| analyzed the planar non-linear a model with random 
dipolar impurities, assuming that the dipole moments are 
annealed and expanding in x /T. In Ref. § we gener- 
alized that analysis, treating the dipoles as quenched. 
The two calculations coincide to lowest order in x, but 
our renormalization group analysis allows a calculation 
of Tn all the way down to zero at x c , supplying a good 
interpolation between these two limits. 

In what follows we summarize that theory, with em- 
phasis on the changes necessary for including the Zn va- 
cancies. We argue that the main effects of the vacancies 
enter in two related ways. First, the concentration z of 
the Zn vacancies renormalizes the concentration of frus- 
trated bonds; when a Cu ion is missing from (at least) one 
end of a "frustrated" bond, then this bond is no longer 
acting like a "dipole" . The probability to find a bond 
without vacancies on both ends is (1 — z) 2 , and therefore 
the effective concentration of "dipolar" bonds is equal to 



y = x(l - zf 



(1) 



Second, when one Cu ion at an end of a hole-doped bond 
is replaced by Zn, then the strong antiferromagnetic cou- 
pling between the spins of the second Cu and of the hole 
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on the oxygen will form a singlet, which is equivalent to a 
magnetic vacancy also on the second Cu. Hence the holes 
increase the number of vacancies, turning their effective 
concentration into 



v = z[l + 2ar(l - z) 



(2) 



In what follows, we shall concentrate on the regime x < 
0.03, where the ^-dependence of v on x is very weak. 

Following Ref . || , we descibe the system by the Hamil- 
tonian 



TL — TL V + TLd , 



(3) 



where H v is the non-linear sigma model (NLctM) Hamil- 
tonian in the renormalized classical region |p^| , represent- 
ing the long wave length fluctuations of the unit vector 
n(r) of antiferromagnetism. In the presence of short- 
range inhomogeneity, this Hamiltonian can be written as 



(4) 



Here i — l,...,d and p = 1, ...,jV run over the spatial 
Cartesian components and over the spin components, re- 
spectively, di = d/dxi, and the effective local stiffness 
p s (r) is a random function. The spatial fluctuations 
Sp s (r) of this function are 5-correlated: [Sp B (r)Sp s (r')] = 
K5(r — r'), where [...] means quenched averages. Simple 
A power counting arguments show Q that K is irrele- 
vant in the renormalization group sense. Therefore, we 
can replace p s (r) in Eq. (^) by its quenched average 
Ps(v) = [p s (r)\. A 

Hd is constructed jl5| to reproduce the dipolar canting 
of the spins at long distances. Denoting by a(r^) the unit 
vector directed along the frustrating bond at r^, and by 
Mim(ri) the corresponding dipolc moment (where m(re) 
is a unit vector giving the direction of the dipole, and Mi 
is its magnitude), we have 



with 



H d = ps(v) J dr^fi(r) • ^n, 



fi(r) = y^8(r- r e )M e ai(r e )m(rt), 



(5) 



(0) 



where the sum runs only over doped bonds which 
frustrate the surrounding (namely have both Cu ions 
present) . 

As argued in Ref. [|| , the renormalization group proce- 
dure generates an effective dipole-dipole interaction be- 
tween the dipole moments, {m(r)}, which is mediated 
via the canted Cu spins. At low temperature T these 
moments develop very long ranged spin-glassy correla- 
tions, and may thus be considered frozen. Hence, we 



treat all the variables r^, a(r^), and m(r^) as quenched, 
and we have 



[/^(r)/ jV (r')] =AVA;<5(r 



(7) 



Here A = Ay, A = M 2 Q/d, where Q = [mi(r)], and the 
effective dipole concentration y replaces the parameter x 
used in Ref. ||. 

With these assumptions, we have now mapped our 
problem to that treated in Ref. ||. We can thus take 
over the results from there, and Tn(x, z) should be equal 
to the Neel temperature derived there for hole concen- 
tration y and stiffness constant p s (v). 

The renormalization group analysis of the Hamiltonian 
(^) H found the two-dimensional antiferromagnetic cor- 
relation length £,2D, as function of the two parameters 
t = T/p s and A. The results contain exact exponential 
factors, which give the leading behavior, and approxi- 
mate prefactors. For doped LCO, the results were given 
for two separate regimes: 



&Va = ^iA°- 8 ex P (g) 



for t < A, and 



eWa = C 2 exp(g[l-(l-^) 3 



(8) 



(9) 



for t > A. Here, a — 3.8A is the lattice constant, and 
the coefficients C\ and C2 may have a weak dependence 
on t and on A. In Ref. Q the data on L^-xSr^CuCM 
were fully described by the constant values C\ = 0.74 
and C 2 = 0.5. 

The three dimensional (3D) Neel temperature was then 
derived from the relation 



c2 



1. 



(10) 



representing the appearance of 3D AFLRO due to the 
weak relative spin anisotropy or the weak relative inter- 
planar exchange coupling, both contained in the param- 
eter a. Combining Eqs. (g) and ([Io|) thus yields an 
a-dependent value for the critical value A c , above which 
AFLRO is lost. This value should give the critical line 
for all t < A. Using the undoped value a ~ 10~ 4 0, Ref. 
Q estimated A c « 0.37. Assuming that a is independent 
of either y and v, this yields y c = X c /A w 0.019, where 
we have used the value A = 20 found for slightly doped 
LCO §. Combining this with Eq. (0), we thus find 



0.019 



(11) 



showing an increase of the antiferromagnetic regime with 
increasing z. At z = 0.15, this would predict x c k, 0.026, 
slightly lower than the observed value. This discrepancy 
could result from various sources. For example, dilution 
may affect the nearest neighbor exchange energy in the 
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plane, J, more strongly than the interplanar interaction 
or the anisotropy. This would imply that a increases 
with z. 

Combining Eqs. (0) and (|l0|), one finds that for t > A, 
the critical line is given by 



BX 



t N (X) „ 

t N (0) " 1 - (1 - 3BX) 1 / 3 



where 



B 



-in 



ln(aC 



tjv(O) 



(12) 



(13) 



We next look at the dependence of Tjy on x for fixed 
z. Ignoring the weak dependence of v on x in Eq. (|^), 
p s is assumed to depend only on z (the relative error 
in p s from neglecting x in Eq. (Q) is less than 3%, see 
below). In that case, p s drops out of the ratio on the LHS 
in Eq. (12), which becomes equal to Tjv(x, z)/Tn(0, z). 



The RHS of that equation now depends only on A = 
Ay = Ax(l — z) 2 , reflecting a universality of the plot of 
Tjv(x, z)/T/v(0, z) versus the rescaled variable y = x(l — 
z) 2 . Note that this universal plot, which should describe 
the Neel temperature for many values of z, requires no 
new parameters; all the parameters are known from the 
limit z = 0. In fact, it is worth noting that the RHS of 
Eq. ( p"2j ) depends only on the combination BX, so that it 
should also apply to other lamellar systems with different 
values of B, resulting from different values of a. 




FIG. 1. The universal plot of Tn(x, z)/Tn(0, z) versus y, 
based on Eqs. (11) (for t < A) and (12) (for t > A). Symbols 
are from experiments: z — 0: empty circles Ref. [4], full 
circles Ref. [19], crosses Ref. [12], triangles Ref. [3]; z — 0.15: 
stars Ref. [4] . 

Figure 1 presents the universal plot of 
T N (x lZ )/T N (0,z), from both Eqs. (O) (for t > A) and 



( |ri| ) (for t < A) . This theoretical curve is then compared 
with various experiments, for both z = and z = 0.15. 
It is satisfactory to note that except for one point, the 
data from the latter are indistinguishable from those for 
the non-dilute case, confirming our universal prediction. 

For comparison of the z— dependence of Tn(x, z) with 
experiments, it is more convenient to scale Tn(x,z) 
by Tn(0,0). For that purpose, we need the ratio 
Tjv(0,z)/XAr(0,0). Theoretically, Eq. @ yields 

T N (0,z) _ t N (0,z)p s (z) _ B(Q) Ps (z) 



7V(0,0) tjv(0,0)p.(0) B(z) Pa (0y 



(14) 



where the weak z— dependence of B{z) may result from 
such a dependence of either a or C2 in Eq. (|l3|). Ac- 
cording to Refs. B and |l7]], the experimental data fit 
the linear dependence 



T N {0,z) 



1 - 3.20z. 



(15) 



^ (0,0) 

up to z — 0.25. At low concentrations, z < 0.10, this 
is in good agreement with the 1/5 expansion result, |l§| ] 
p s (z)/p s (0) = 1 — 3.14z, if one uses the approximation 
B(z) « B(0) in Eq. (pT[). At higher concentrations the 
ratio p s (z) / p s (0) in the classical limit decreases with di- 
lution approximately as p s (z)/p s (0) = 1 — 3.14z + 1.57z 2 , 
]l8| i.e. slower than the experimental Tn{z)/T^{Q). This 
discrepancy can be due to quantum corrections to p s , or 
to the z dependence of B. Substituting Eqs. (|l4| ) and 
( |l5| ) (also replacing z by v) into Eq. ([lj), we have 



T N (x,z) 
TV (0,0) 



(1 - 3.20w) 



ABy 



1 - (1 -SABy) 1 /*' 



(16) 
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FIG. 2. T N (x, z)/T N (0,0) versus z. The full, dashed, and 
dash-dotted line represent Eq. (16), for x — 0.018, x = 0.017, 
and x = 0.016 respectively. Circles are experimental data, 
from Ref. [4], for Sr concentration 0.017 ± 0.001. 
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Figure 2 shows the dependence of Tjv(x, z) on the di- 
lution z, given by the above equation, for three concen- 
trations of holes. The theoretical curve for x = 0.018 
reproduces very well the observed maximum in the de- 
pendence of Tn(x, z)/Tn(0, 0) on z. The experimen- 
tal points were measured at a nominal Sr concentration 
x = 0.017 ±0.001. 

An important prediction of the theory is the high sen- 
sitivity of the maximum to the hole concentration. The 
maximum exists only at x sufficiently close to x c , and 
disappears at lower x. It would be interesting to check 
this prediction experimentally. 

In conclusion, we found the combined effect of hole 
doping and magnetic dilution on the long-range order in 
lamellar Heisenberg antiferromagnets. We showed that 
dilution weakens the destructive effect of the holes on 
the AFLRO. The critical concentration x c increases with 
dilution, and the dependence of Tjv on vacancy concen- 
tration reveals a maximum, if the hole concentration is 
sufficiently close to x c . These findings are in quantita- 
tive agreement with the experiment. Furthermore, the 
experimental data for the two concentrations agree with 
each other even better than with the theoretical curve, 
demonstrating the validity of the scaling x — > x(l — z) 2 
beyond any theory. We note that the consistency of our 
theory with the experiments, with no new adjustable pa- 
rameters, supports the validity of the frustration model. 
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there for the NSF under grant No. PHY94007194. 
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